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ABSTRACT 

Power density spectra (PDSs) of long gamma-ray bursts (GRBs) provide useful information on GRBs, in- 
dicating their self-similar temporal structure. The best power-law PDSs are displayed by the longest bursts 
(Tgo > 100 s) in which the range of self-similar time scales covers more than 2 decades. Shorter bursts have 
apparent PDS slopes more strongly affected by statistical fluctuations. The underlying power law can then be 
reproduced with high accuracy by averaging the PDSs for a large sample of bursts. This power law has a slope 
ol ~ —5/3 and a sharp break at ~ 1 Hz. 

The power-law PDS provides a new sensitive tool for studies of GRBs. In particular, we calculate the PDSs of 
bright bursts in separate energy channels. The PDS flattens in the hard channel (hv > 300 keV) and steepens in 
the soft channel (hv < 50 keV), while the PDS of bolometric light curves approximately follows the —5/3 law. 

We then study dim bursts and compare them to the bright ones. We find a strong correlation between the burst 
brightness and the PDS slope. This correlation shows that the bursts are far from being standard candles and 
dim bursts should be intrinsically weak. The time dilation of dim bursts is probably related to physical processes 
occurring in the burst rather than to a cosmological redshift. 

Subject headings: gamma rays: bursts 



1. INTRODUCTION 



The results are discussed in 57. 



The light curves of gamma-ray bursts (GRBs) typically have 
many random peaks and in spite of extensive statistical studies 
(e.g., Nemiroff et al. 1994; Norris et al. 1996; Stern 1996), the 
temporal behavior of GRBs remains a puzzle. Contrary to the 
complicated diverse behavior in the time domain, long GRBs 
show a simple behavior in the Fourier domain (Beloborodov, 
Stern, & Svensson 1998). Their PDS is a power law of in- 
dex a w —5/3 (with a break at ~ 1 Hz) plus standard (expo- 
nentially distributed) statistical fluctuations superimposed onto 
the power law. The PDS slope and the break characterize the 
process randomly generating the diverse light curves of GRBs. 
Intriguingly, the PDS slope coincides with the slope of the Kol- 
mogorov law. The power-law behavior is seen in individual 
bursts (we illustrate this in §3) and may provide a clue to the 
nature of GRBs. 

In the present paper, we study in detail the PDSs of gamma- 
ray bursts. In our analysis we use 527 GRB light curves with 
64 ms resolution obtained by the Burst and Transient Source 
Experiment (BATSE) in the four LAD energy channels, I-IV: 
(I) 20 - 50 keV, (II) 50 - 100 keV, (III) 100 - 300 keV, and (IV) 
hv > 300 keV. The background is subtracted in each channel 
using linear fits to the 1024 ms data. 

The method of data analysis is described in §2. In §3, we 
study a sample of the 4 brightest and longest bursts. In §4 the 
average PDS for the full sample of 527 GRBs is calculated and 
discussed. In §5, we calculate the PDSs in the separate energy 
channels and quantify the difference of the temporal structure 
between the channels in terms of the PDS slope. We also calcu- 
late the autocorrelation function (ACF) in the separate channels 
and compare the results with previous studies of the ACF. In 
§6, we address dim GRBs and compare them to the bright ones. 



2. DATA ANALYSIS 

2.1. The Sample 

Long bursts are of particular interest since their internal 
temporal structure can be studied by spectral analysis over a 
larger range of time scales. We choose bursts with durations 
Tgo > 20 s where Tgo is the time it takes to accumulate from 
5% to 95% of the total fluence of a burst summed over all the 
four channels, I+II+III+IV. Hereafter, we measure the bright- 
ness of a burst by its peak count rate, C pea k, in channels II+III. 
We analyze bursts with C pC ak > 100 counts per time-bin, 
At = 0.064 s. This condition coupled with the duration condi- 
tion, Tgo > 20 s, gives a sample of 559 GRBs. The sample still 
contains GRBs with low fluence, which are not good for Fourier 
analysis. We exclude bursts with fluences $ < 32C pC ak- The 
resulting sample contains 527 bursts. 

2.2. PDS Calculation 

We calculate the Fourier transform, C/, of each light curve, 
C(t), using the standard Fast Fourier Transform method. The 
power density spectrum, Pf, is given by Pf = (CfCj + 
C_/C* / )/2 = C S C) since C{t) is real. The Fourier trans- 
form is calculated on a standard grid with a time bin At = 64 
ms and a total number of bins Nbm = 2 14 , which corresponds 
to a total time T w 1048 s since the trigger time. The light 
curve of each burst is considered in its individual time window 
(t\,t2) (see §2.4). In the calculations of the PDS, we extend 
the time interval to (0, T) by adding zeros at (0, t\) and (t2,T). 
The adding zeros introduces random small-scale fluctuations in 
the PDS (see, e.g., Bracewell 1967). We, however, study the 



PDSs averaged over adjacent frequencies and/or over a sample 
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FIG. 1. — The light curves (in channels II+III) and their PDSs for the four brightest bursts in the sample with T90 > 100 s. The light curves are taken peak- 
normalized and, correspondingly, the PDSs are normalized by Cp Cak . Since the bursts are very bright, their Poisson levels (normalized by Cp Cak ) are very low, 

Po ~ 10~ 4 . The straight lines show the fit to the average of the 4 PDSs, log Pf = 1.25 — 1.45 log / (see Fig. 2). 
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of GRBs. Then the fluctuations in Pf associated with a specific 
choice of the grid disappear and do not affect the results. 

2.3. The Poisson Level 

Poisson noise in the measured count rate affects the light 
curve on short time scales, and, correspondingly, affects the 
PDS at high frequencies. The Poisson noise has a flat spectrum, 
Pf oc f°, introducing the "Poisson level", P , in a PDS. This 
level equals the burst total fluence including the background in 
the considered time window. The power spectrum above this 
level displays the intrinsic variability of the signal (but see §4.2 
for the time-window effects). We calculate the individual Pois- 
son level for each burst and subtract it from the burst PDS. 

2.4. The Time Window 

One would like to see the whole burst in the window. How- 
ever, (1) it is difficult to determine exactly the end of a burst 
because the burst can always have a weak "tail" hidden in the 
background; (2) the window should not be very large because 
inclusion of long weak tails leads to an increase of the back- 
ground fluence, without a substantial increase in the signal 
fluence, $. The resulting low ratio implies a high Poisson 
level, which makes the quality of the PDS worse. 

To reduce the Poisson level, we cut off the light curves at a 
time i 2 defined so that the signal count rate C(t) does not ex- 
ceed eCpeak at t > t 2 and C(t 2 ) — eC pea k- Keeping in mind 
bursts with a weak beginning, we define the starting window 
time, ti, so that C(t) < eC pea k at t < t\ and C(ti) = £C pca k- 
In our analysis, e = 0.05 is chosen. We checked that vary- 
ing e does not significantly affect the results unless e > 0.1. 
Note that, for dim bursts, Poisson fluctuations of the back- 
ground (which are imprinted on C (t) even after subtraction of 
the average background level) may exceed £C poa k. Then the 
window (ti, i 2 ) is determined by the background fluctuations 
rather than by the signal. 

3. INDIVIDUAL BRIGHT BURSTS 

The brightest and longest bursts are the best ones for Fourier 
analysis. In this section, we study the four brightest bursts 
with T 90 > 100 s. They have trigger numbers # 2156 (GRB 
930201), 2856 (GRB 940302), 3227 (GRB 941008), and 6472 
(GRB 941110). 

We take the light curves, C(t), summed over channels II and 
III, in which the signal is strongest. To simplify the comparison 
of different bursts, we take peak-normalized light curves. Their 
Fourier transform, C/, is therefore normalized by C pca k, and 
the PDS, P f = C f C), is normalized by C£ eak . 

The four peak-normalized light curves and their PDSs are 
shown in Figure 1 (we smooth the PDSs on the scale A log / = 
0.04 before plotting). The light curves are very different while 
the PDSs are similar. They can be described as a single power 
law, logP/ = A + a log / (Awl and a s=y —1.5) with super- 
imposed fluctuations APf/Pf <~ 1. In spite of the large APf, 
the power-law behavior can be seen in each burst due to the 
power-law extending over more than two decades in frequency. 

The presence of an underlying power law is an interesting 
feature of the GRB temporal behavior. A possible way to ex- 
tract it from the noisy individual PDSs is to take the average 
PDS over a sample of long GRBs. Then the fluctuations af- 
fecting each individual PDS tend to cancel each other and one 
can see the power-law behavior. The average of the four PDSs 
is plotted in Figure 2. One can see that the amplitude of ir- 
regularities is reduced after the averaging and one can measure 



the slope of the resulting PDS. The power-law fit in the range 
-2.0 < log / < 0.2 gives log P f = 1.25 - 1.45 log/. (We use 
the standard \ 2 fitting procedure with equal weights for each 
frequency bin.) 
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FIG. 2. — The average of the 4 PDSs shown in Figure 1. The line shows the 
power-law fit, log Pf = 1.25- 1.45 log /, in the range -2.0 < log / < 0.2. 

The PDS analysis was previously performed for a number of 
individual bursts. Belli (1992), for example, analyzed GRBs 
detected by the Konus experiment onboard the Soviet probes 
Venera 1 1 and Venera 12, and Giblin et al. (1998) studied indi- 
vidual BATSE bursts. The power-law PDS can be observed in 
their longest complex bursts as well. 

When turning to GRBs with relatively modest durations, 
T90 ~ 20 — 100 s, it becomes more difficult to see the power- 
law behavior in individual bursts. The statistical fluctuations, 
APf/Pf ~ 1, significantly affect the apparent PDS slope de- 
rived for the available range of frequencies, which is limited by 
<~ T 9 q . However, having a large number of bursts allows one 
to study the PDS by averaging over the sample. 

4. THE AVERAGE PDS 

In most of the GRB models (e.g., in the internal shock 
model), different bursts are produced by one physical mecha- 
nism of stochastic nature, i.e., an individual burst is a random 
realization of the same standard process. Stern & Svensson 
(1996) suggested an empirical pulse-avalanche model of this 
type. They showed that near-critical avalanches reproduce the 
statistical characteristics of GRBs as well as their extreme di- 
versity. With such an approach, individual GRBs are like pieces 
of a "puzzle", and the features of the standard engine can be 
probed with statistical methods applied to a sufficiently large 
ensemble of GRBs. 

4. 1 . The averaging 

The simplest statistical characteristic of the PDS is the aver- 
age PDS, Pf. The averaging means that we sum up the PDSs 
of individual bursts with some weights (normalization) and di- 
vide the result by the number of bursts in the sample. If no 
normalization is employed, then the brightest bursts strongly 
dominate Pf, and their individual fluctuations lead to strong 
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fluctuations in Pf. The PDSs of relatively weak bursts are sup- 
pressed proportionally to C~ e ^ k and they are lost in the averag- 
ing. A normalization procedure is thus needed to increase the 
weight of relatively weak bursts in the sample. We prefer the 
peak-normalization because: 

(i) The amplitude of the resulting PDS is practically indepen- 
dent of the burst brightness (see, e.g., Fig. 1). 

(ii) The fluctuations in Pf turn out to be minimal and the best 
accuracy of the slope is achieved. The distribution of the indi- 
vidual Pf around Pf follows a standard exponential law (see 
Beloborodov et al. 1998), and the amplitude of fluctuations in 
Pf is given by a simple formula, APf/Pf <~ iV -1 / 2 , where N 
is the number of bursts in the sample. 

The averaging procedure is the following. (1) For each burst 
we first determine its individual Poisson level. (2) We deter- 
mine the peak, C pca io of the light curve, C(t), with the back- 
ground subtracted. (3) We normalize the light curve to its peak. 
(4) We calculate the PDS of the normalized light curve. (5) We 
normalize the Poisson level by Cp eak and subtract it from the 
PDS. (6) We sum up the resulting PDSs over the sample and 
divide by the number of bursts. The resulting average PDS, Pf, 
for the 527 peak-normalized light curves in channels (II+III) is 
shown in Figure 3. 

For comparison, we also plot the average PDS with the 
normalization, where $ is the burst fluence excluding the back- 
ground. Each light curve is then normalized by ^/W, and its 
PDS (and its Poisson level) is normalized by $. This normal- 
ization gives different weights (the weights of dim bursts are 
reduced as compared to the peak-normalization). The resulting 
average PDS has a different amplitude. Nevertheless, Pf again 
follows a power law with approximately the same slope. This 
provides evidence that the self-similar behavior with a m —5/3 
is an intrinsic property of GRBs, rather than an artifact of the 
averaging procedure. This interpretation is also supported by 
the fact that we observe the power law in the longest individual 
bursts. When comparing individual PDSs, Pf, with the peak- 
normalized Pf, one sees that Pf is exponentially distributed 
around Pf, and the distribution is self-similar with respect to 
shifts in /. 

The power-law fit to Pf in the range —1.4 < log/ < 
—0.1 gives a s=a —1.75 for the peak-normalized bursts, and 
a « —1.67 for \/$-normalized bursts. The formal error found 
from the \ 2 fitting is very small, Aa <~ ±0.01 (we take the 
variance, APf/Pf = N~ 1/2 ). Note that the slope in the 
peak-normalized case is different from —1.67 reported in Be- 
loborodov et al. (1998). This change is caused by the fact that 
we have extended the sample to smaller brightnesses (see §6). 

The deviation from the power law at the low-frequency end 
is due to the finite duration of bursts. At the high frequency 
end, there is a break at <~ 1 Hz. The break is observed in the 
brightest GRBs even without subtracting the Poisson level (see 
Beloborodov et al. 1998 and the top panel in Fig. 8). Note that 
the break position is the same for the peak-normalized and \/$- 
normalized bursts (see Fig. 3b). It stays the same in the separate 
energy channels (see Fig. 5) and does not depend on the Pois- 
son level. One may also see the break in individual long bursts 
(Figs. 1 and 2). The break is far too sharp to be explained as an 
artifact of the 64 ms time binning, which suppresses the PDS 
by a factor of [sin(7r ■/ At) / '(tt / 'At)} 2 where At = 64 ms is the 
time bin (cf. van der Klis 1989). 
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FIG. 3. — The average PDS (in channels II+III) for the full sample 
of 527 GRBs. (a) - The histograms marked as 1 and 2 correspond to the 
peak-normalization and V$-normalization, respectively. The horizontal lines 
show the corresponding normalized Poisson levels averaged over the sample. 
Dashed lines are the power-law fits in the range —1.4 < log/ < —0.1: 
log Pf = 0.96 — 1.75 log / in the peak-normalized case and log Pf = 
1.76 — 1.67 log/ in the \/i>-normalized case. Dotted line shows the —5/3 
slope, (b) - Same as (a) except that now Pf and the Poisson levels are plotted 
multiplied by / 5 / 3 . 



4.2. The Effects of Finite Signal Duration 

Fourier analysis was designed for physical problems dealing 
with linear differential equations. For example, it is usually ap- 
plied to small perturbations above a given background solution. 
The Fourier power spectrum is also commonly used in the tem- 
poral studies of long signals or noises, e.g., in persistent astro- 
physical sources. By contrast, GRBs have strongly non-linear 
signals with short durations. The typical number of BATSE 
time-bins (At = 64 ms) in a long GRB is a few x 10 3 . Do the 
effects of finite duration (i.e., time-window effects) strongly af- 
fect the measured PDS? 
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FIG. 4. — The average PDS, Pf, for the 64 fragments of noise with a power- 
law PDS of index a (see text for details). Dotted lines show the original PDS 
slope. The deviation of Pf from this slope is due to the time-window effects. 

The issue is illustrated in Figure 4. We prepared an artificial 
long signal with Pf oc /" and exponentially distributed Pf, 
i.e., the probability to detect Pf at a given / is proportional to 
exp(— Pf /Pf). The phase structure was taken to be Gaussian 
(i.e., random). The signal duration is T = 2 17 Ai w 8400 s. 
Then, we cut the long signal into 64 pieces of equal length 
to = 2 n Ai w 130 s. We thus get 64 short signals, each is 
a random realization/fragment of the same stationary process 
characterized by the index a. Analogously to our analysis of 
real GRBs, we normalize each signal to its peak and add ze- 
ros up to T — 2 14 Ai (our standard grid). Then we calculate 
the average PDS, Pf, for the 64 artificial signals. The result is 
shown in Figure 4 for the three cases: a = (Poisson), a = — 1 
(flicker), and a — —5/3 (Kolmogorov). 

One can see that the time-window effects are strong in the 
case a = 0. The Poisson signal is roughly constant on large 
time scales. As a result, at modest /, Pf is just the power spec- 
trum of the to-window. By contrast, in the case a = —5/3, 
we have large-amplitude variations in the signal on all time 
scales. The average PDS of the 64 short signals then repro- 
duces well the intrinsic PDS slope throughout the whole range 
of frequencies, down to / ~ t$ (see bottom panel in Fig. 4). 
Hence, the power law we observe in the average PDS of GRBs 
can be interpreted as that GRBs are random short realizations 
of a standard process which is characterized by the PDS slope 
a w —5/3. Note, however, that the power spectrum does not 
provide a complete description of the signal since the phase 
structure is not considered. 



152 bursts. 

The burst analysis is now performed in each channel sepa- 
rately. We determine the Poisson level of a burst in each chan- 
nel, Pq, and find the peak of the light curve, Cp Cak , where 
i = I, II, III, IV. Then we perform the peak-normalization: 

C*(t) C*(i)/^ oak and P* ^ P 7(^ cak ) 2 . We set the 
time window in each channel (t\, t\) as described in §2.4. 




FIG. 5. — The average PDS for the sample of 152 GRBs in separate energy 
channels, I, II, III, and IV. Dashed lines show the power-law fits (see Table 1). 
Solid lines show the Poisson level. 



5. THE PDS AND THE ACF IN CHANNELS I, II, III, IV 

What does the average PDS look like in the separate en- 
ergy channels? The signal to noise ratio is low in channel I 
and especially in channel IV. The number of GRBs for which 
good PDSs can be obtained in channels I-IV is therefore lim- 
ited to the brightest bursts. We choose a sample of bursts with 
Cpeak > 500 counts/bin (in channels II+III), which contains 



The resulting average PDSs in channels I-IV are shown in 
Figure 5. The differences in the slopes are clearly seen. We 
fitted the PDSs by power laws, logP/ = A + a log /, in the 
range —1.6 < log/ < 0. Channel IV is fitted in the range 
-1.3 < log / < -0.1. A and a of the fits are listed in Table 1. 
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Table 1. 



Channel 

a 
A 



I 

-1.72 
1.03 



II 

-1.67 
1.05 



III IV 

-1.60 -1.50 
1.06 1.07 



The observed behavior can be briefly described as follows: 
The "red" power (Pf at low /) decreases at high photon ener- 
gies and increases at low photon energies. It should be com- 
pared with the well-known fact that the pulses in a GRB are 
narrower in the hard channels (e.g., Norris et al. 1996). The 
hardness of emission varies dramatically during a burst and this 
leads to different temporal structure in different channels. E.g., 
pulses observed in the soft channels may be suppressed in the 
hard channels. One therefore could expect that the PDS has 
different slopes in different channels. Note that, typically, most 
of the GRB energy is released in channels II+III, and the aver- 
age PDS of bolometric light-curves approximately follows the 
-5/3 law. 

In principle, the autocorrelation function (ACF) contains the 
same information as the PDS, since one is the Fourier transform 
of the other (the Wiener- Khinchin theorem). In practice, the 
two are not completely equivalent because of the time-window 
effects and the presence of a noisy background. On modest 
time scales, 30 s, the direct ACF calculation and the calcula- 
tion via the Fourier transform of the PDS give the same result 
to within a few percent (we have calculated the ACF by both 
methods). The average ACF, A(t), for our sample of 152 bright 
GRBs is shown in Figure 6 for each of the four channels. The 
ACF gets narrower at high energies, in agreement with previous 
studies (Fenimore et al. 1995), except that our ACF is <~ 2 times 
wider as compared to that obtained by Fenimore et al. (1995). 
The ACF width averaged over the channels is in approximate 
agreement with that calculated for the bolometric light curves 
by Stern & Svensson (1996). 

Note that the average ACF is obtained by summing up the 
individual ACFs normalized to unity, i.e., A(0) — 1. Recalling 
that the ACF is the Fourier transform of the PDS, one can see 
that this normalization is equivalent to the normalization of the 
light curve by V^U where P tot = J[C(t)] 2 dt = J P f df is 
the total power. This normalization is different from the peak- 
normalization we use in the calculations of the average PDS, 
i.e., we prescribe different weights to individual bursts when 
averaging the PDS and the ACF, respectively. Therefore, the 
average ACF is not the Fourier transform of the average PDS 
shown in Figure 5. This relation holds when the average PDS 
is calculated with the V-Ptot -normalization. 

The ACF in each channel is perfectly fitted by the stretched 
exponential: A(t) = cxp(— [t/t ] /3 ) (see Fig. 7). The parame- 
ters ro and (3 are listed in Table 2. The index (3 is related to the 
PDS slope a by the simple relation (3 ~ — (1 + a). 



Table 2. 



Channel 

TO 



I 

0.73 
14.0 



II 

0.67 
10.7 



III IV 

0.63 0.6 
7.3 5.1 



The changing (3 demonstrates that the ACF shape changes 
from channel to channel, as it should do since the PDS changes. 



As a first approximation, the PDS slope is equal to —5/3, and 
the ACF index (3 « 2/3 = 5/3 — 1. 




o 
<c 

w 
o 



-0.5 - 



, ' ' ' ' 


1 1 1 1 1 1 


- 




















\ \ ^ 






\ x X 

\ \ X 
\ \ 


I 




\ \ 


- 
















X, \ 






iv X \^ 










II 


i i i i i i i i 


i 



10 



20 



30 



T, S 



FIG. 6. — The average autocorrelation function for the sample of 152 bright 
bursts, in channels I-IV. 
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FIG. 7. — The average ACFs in channels I-IV are plotted against t 13 which 
should give a straight line for a stretched exponential. 

The values of to quantify the ACF width in different chan- 
nels. The parameter r is better defined as compared to mea- 
suring the ACF width at a certain level, e.g., at e~° 5 of the 
maximum level, as done in Fenimore et al. (1995). Neverthe- 
less, the scaling of tq with energy is approximately the same, 
To oc E~° A where E is the photon energy of the low energy 
boundary of the channel, see Fenimore et al. (1995). It should 
be stressed, however, that the stretching of t is not related to 
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the stretching of any physical time scale of intrinsic correla- 
tions during the burst. The ACF width, t , is rather related 
to the position of the breaks in the PDS, especially the low fre- 
quency break, which in turn is determined by the burst duration. 
For an infinitely long signal with a power law PDS, To would 
tend to infinity. This is natural since the power law behavior 
implies self-similarity, i.e., the absence of any preferred time 
scale. Specific time scales are introduced only by the breaks 
in the PDS. We therefore have only two physical time scales in 
long GRBs: The first one is associated with the 1 Hz break and 
the second one is associated with the low frequency break due 
to the finite burst duration. 

6. DIM VERSUS BRIGHT BURSTS 

Now we address the following question: Is there any corre- 
lation between the PDS slope and the burst brightness? 



6.1. PDS Slope Correlates with the Burst Brightness 

We take the full sample of 527 light curves in channels 
II+III and divide it into 3 groups of different brightnesses: (A) 
Cpeak > 800, (B) 300 < C pcak < 800, and (C) 100 < C poak < 
300. Group A contains 91 bursts, group B 222 bursts, and group 
C 214 bursts. We have calculated the average PDS for each 
group separately employing the peak-normalization. The re- 
sults are presented in Figure 8. We fitted the average PDSs by 
power laws, log P/ = A + a log /. The parameters A and a of 
the fits in the range — 1.5 < log / < —0.1 are listed in Table 3. 
We conclude that the average PDS of dim bursts gets markedly 
steeper. 

Table 3. 



group ABC 

a -1.63 -1.74 -1.82 
A 1.04 0.98 0.91 



6.2. Subtraction of the Poisson Level 

In the PDS calculations we subtracted the Poisson level, 
which is quite high for dim bursts (see Fig. 8). One therefore 
should address a technical question: How well is the intrinsic 
PDS restored after subtracting the Poisson level? 

To investigate this issue we have performed the following 
test. We take the sample of 91 bright bursts (group A) and 
add a strong Poisson noise (with an average level of 25000 
counts/bin) to each burst in the group. We thus artificially 
increase the Poisson level by two orders of magnitude. Al- 
ternatively, one may consider this procedure as a rescaling of 
bright bursts to a smaller brightness while keeping the Poisson 
background at the same level. Then, we analyze the artificial 
bursts in the same way as we did with the real dim GRBs. Note 
that even after subtraction of the time-averaged background, its 
Poisson noise strongly affects the signal and may create artifi- 
cial peaks dominating the true peak of the signal. For the arti- 
ficial bursts, we know the true peak (which is the peak of the 
original bright burst). In real dim bursts we do not know the 
position of the true peak and employ the peak search scheme 
described in Stern, Poutanen, & Svensson (1999). 
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FIG. 8. — The average PDSs for the three brightness groups A, B, and 
C. Dashed lines show the power-law fits (see Table 3). Solid lines show the 
Poisson level. 

The result is compared with the average PDS of the original 
sample in Figure 9. We find that the subtraction of the Poisson 
level allows one to restore the original PDS, P/, also at fre- 
quencies where Pf is well below the Poisson level. Even the 1 
Hz break remains present. We conclude that the high Poisson 
level is unlikely to significantly affect the measured PDS slope. 

7. DISCUSSION 

7.1. Relation to the Average Time Profile 

The average time profile (ATP) of GRBs was found to follow 
a stretched exponential of index 1/3 (Stern 1996). Is there any 
relation between the ATP and the average PDS? One should 
note two important differences between the ATP and the PDS 
studies: (1) the —5/3 PDS, though affected by the statistical 
fluctuations, is observed in individual bursts, while the ATP is 
a purely statistical property of a large sample of bursts; (2) the 
ATP is constructed for GRBs of all durations (and only in this 
case it displays the perfect stretched exponential) while the PDS 
is studied for long bursts only. Nevertheless, both the ATP and 
the PDS characterize the stochastic process generating GRBs. 
Stern (1999) shows that both can be reproduced simultaneously 
with the pulse-avalanche model of Stern & Svensson (1996). 
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FIG. 9. — The average PDS of 91 artificial dim bursts created from group 
A of the brightest bursts (see Fig. 8). Solid histogram - the original PDS for 
group A. Dashed and dotted histograms - the average PDS of the artificial 
bursts with and without subtraction of the Poisson level, respectively. The two 
solid lines show the average Poisson levels for the original and the artificial 
bursts, respectively. 



7.2. The 1 Hz Break 

The sharpness of the break in the average PDS appears to be 
an important feature that constrains models of GRBs. If the sig- 
nal is produced in the rest frame of a relativistic outflow then 
variations of the outflow Lorentz factor, Y, from burst to burst 
would smear out the break. The sharpness of the break then 
implies a narrow dispersion of V, AY /Y 2, which appears to 
be unlikely. 

Alternatively, the GRB variability may come from the central 
engine. The break in the PDS might then be related to a typical 
time scale, <~ 1 s, in the central engine. 

Finally, one may associate the break with the dynamical time 
scale corresponding to the inner radius of the optically thin zone 
of the outflow, Then the variability on time scales shorter 
than t* = R*/cY 2 can be suppressed. In this case, however, 
one should explain why t* stays the same in different bursts. 

7.3. Dim GRBs 

It has often been hypothesized that dim bursts are at high 
cosmological redshifts. For instance, it is necessarily the case 
if GRBs have approximately the same intrinsic luminosity, the 
so-called "standard candle" hypothesis. We can test this hy- 
pothesis using the power spectrum analysis. 

Suppose that the dim bursts are intrinsically the same as the 
bright ones. Then any difference in their average PDS should 
be due to a cosmological redshift. First consider the bolometric 
light curves assuming that their average PDS follows the —5/3 
power law. It is easy to see that a redshift, z, will not change 
the PDS slope. As we normalize each bursts to its peak before 
the averaging, the effect of a redshift is just stretching the light 
curve in time, i.e., precisely the time dilation effect. This will 
lead to an increase of the net normalization of the PDS by a fac- 
tor of (1 + z) 1 / 3 . The slope does not change since the dilation 
factor (1 + z) is the same for each time scale. 



One should, however, recall that we observe bursts in a lim- 
ited spectral interval. A redshift then implies a shift of the sig- 
nal with respect to our spectral window. E.g., photons detected 
in channel III would originally have been emitted in channel 
IV. As we know from §5, the PDSs are different in different en- 
ergy channels. Therefore, one expects that the PDS slope will 
change for redshifted GRBs. 

We have seen in §5 that the PDS of bright bursts flattens in 
the hard channels. Hence, a redshift of the bright bursts must 
be accompanied by a flattening of the PDS. Contrary to this be- 
havior, we observe that the PDS of dim bursts steepens. Hence, 
the evolution of the PDS with brightness is inconsistent with the 
standard candle hypothesis. It implies that the burst luminosity 
function is broad and dim bursts are intrinsically weak. 

Evidence for a broad luminosity function is also found when 
looking at the isotropic luminosities of the bursts with measured 
redshifts. Note, however, that the differences in the apparent 
luminosities could be caused by orientation effects if GRBs are 
beamed. One should not therefore exclude that the total intrin- 
sic luminosities are approximately the same. The different tem- 
poral structure of dim bursts may be an important fact in this 
respect. In particular, it suggests that the observed time dilation 
of dim bursts (e.g., Norris et al. 1994) may be caused mainly 
by physical processes occurring in the bursts rather than by a 
cosmological redshift. Note that the intrinsic difference of the 
temporal structure of dim GRBs was also found when study- 
ing their average time profile (see Stern et at. 1997). The ris- 
ing part of the ATP does not change with decreasing brightness 
while the decaying part suffers from time dilation. This be- 
havior is inconsistent with a cosmological time dilation which 
should apply equally to both parts of the ATP. 

7.4. Internal Shocks 

A likely scenario of GRBs involves internal shocks in a rela- 
tivistic outflow with a Lorentz factor Y ~ 10 2 (see, e.g., Piran 
1999 for a review). The shock develops when an inner faster 
shell of the outflow tries to overtake the previous slower shell. 
The pulses in a burst are then associated with collisions between 
the shells. The PDSs predicted by the model were recently 
tested against the observed —5/3 law (Panaitescu, Spada, & 
Meszaros 1999; Spada, Panaitescu, & Meszaros 2000; Be- 
loborodov 1999,2000). The self-similar temporal structure with 
a = —5/3 can be reproduced by the model. Further constraints 
should, however, be imposed by the observed dependence of a 
on photon energy. Besides, the phase structure of the signal, ne- 
glected so far, should be taken into consideration (Beloborodov 
2000). 
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